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Quasiclassical approach to the weak levitation of extended states in the quantum Hall 

effect 
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The two-dimensional motion of a charged particle in a random potential and a transverse magnetic 
field is believed to be delocalized only at discrete energies Em- In strong fields there is a small 
positive deviation of from the center of the Ai'th Landau level, which is referred to as the "weak 
levitation" of the extended state. I calculate the size of the weak levitation effect for the case of 
a smooth random potential re-deriving earlier results of Ifaldane and Yang [Phys. Rev. Lett. 78, 
298 (1997)] and extending their approach to lower magnetic fields. I find that as the magnetic 
field decreases, this effect remains weak down to the lowest field Bmin where such a quasiclassical 
approach is still justified. Moreover, in the immediate vicinity of Bmin the weak levitation becomes 
additionally suppressed. This indicates that the "strong levitation" expected at yet even lower 
magnetic fields must be of a completely different origin. 



The levitation of the extended states in the auan- 
tum HalLeffect has been proposed by KhmclnitskiiB and 
Laughlino (KL) almost fifteen years ago. Based on scal- 
ing arguments (see Ref. ^ for review) KL further sug- 
gested that the energies En of the extended states can 
be obtained by solving the equation 



a^yiEp^EN) - (^ + 1/2) 



(1) 



where axy is the unrenormalized (short length scale) 
Hall conductance, Ep is the Fermi energy, and N runs 
through the set of integer values iV = 0, 1, . . .. Using the 
Drude-Lorentz formula 
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for the left-hand side, they obtainedQ 

En = iN + l/2)hu;c [1 + {lu,t)-^] 



(2) 



(3) 



where coc is the cyclotron frequency, r is the zero-field 
momentum relaxation time, and n is the electron den- 
sity. Thus, in strong fields where lUcT 3> 1 the energy 
En of the A''th extended state is close to the center of 
the Nth Landau level E"^ = fiUc{N -f 1/2). As the mag- 
netic field decreases. En floats upward with respect to 



i?^, so that the relative deviation 5E 
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mcreases 



SEn/E^ = {En - E^)/E^ = (u^t)- 



(4) 



I will call the regimes SEn/E"^ < 1 and 5En/E^ > 1 
a "weak" and a "strong" levitation regimes, respectively. 
The regime of strong levitation is, of course, the most 
interesting. Unfortunately, this regime is also the hard- 
est one to study. Up to date there is no progress in 
analytical treatment of this problem; as for the original 
arguments,!!! they suffer from the absencer-pf completely 
satisfactory derivation of the scaling laws.B 



Not long ago, the study of the weak levitation ph&. 
nomenon has been pioneered by ShahbazyanjAnd Raikhc! 
and then continued by several other groups .UB The ad- 
vances in the analytical treatment of the problem are 
due to the existence of a transparent physical picture 
of localization in sufficientlji. strong magnetic field and a 
smooth random potential.^ I will start with briefly re- 
calling this picture, which leads to an approximate equal- 
ity En ~ . Then I will review the arguments of 
Haldane and Yangfl (HY) who identified the leading con- 
tribution to 5En = En — E^. Finally, I will present my 
own results. 

The first simplification of the strong magnetic field 
limit comes from the approximate separation of the elec- 
tron's motion into a fast rotation along the cyclotron 
orbit and a much slower dynamics of its guiding center. 
Unfortunately, the guiding center coordinates (unlike the 
coordinates of the particle itself) can not have definite 
values simultaneously. The characteristic uncertainty in 
the guiding center position is of the order of the mag- 
netic length I = ^h/mujc- And here the smoothness of 
the random potential brings the second crucial simplifica- 
tion: as long as the correlation length d of the potential is 
much larger than I, this uncertainty can be ignored. As a 
result, there are two approximate integrals of motion: the 
energy of the cyclotron motion, i.e., the kinetic energy, 
and the energy of the guiding center degree of freedom, 
which is essentially the potential energy. In this approx- 
imation the guiding center is permanently bound to a 
certain level line U{x,y) — const of the random potential 
L/(x, y) and performs a slow drift along such a line.|-The 
extended states correspond to unbounded level lines.!l2! In 
fact, for a wide class of potentials statistically symmetric 
under the sign change, the unbounded level line_£perco- 
lation contour) is unique and is at zero energy.!li! Since 
the total energy of the electron is equal to the sum of the 
energies of the cyclotron and the drift degrees of freedom. 
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and the cyclotron energy is quantized in hujc quanta, in 
this approximation is equal to . 

The calculation of the the weak levitation correction 
requires taking into account the Landau level mixing. 
HY demonstrated that such a mixing simply modifies the 
form of the potential in which the guiding center drifts. 
The new potential is not statistically symmetric under 
sign change. In fact, its percolation level is at higher 
energy.Q Correspondingly, Ef^ is larger than E"^ . 

Using the quantum-mechanical perturbation theory, 
HY found that 



5En - (iV + l/2){W'^/huJc){l/df 



(5) 



where W is the rms amplitude of the random potential. 
I retained only the first term in the perturbation series 
obtained by HY. The next term is smaller if d ^ IN^^^. 
In addition, HY require that huic ^ W. Denote by Rc 
the classical cyclotron radius at energy E^^f. It is easy to 
see that Rc = \/2N +11, so that Eq. (||) becomes 



6En/E^ ^{W/EN)HRc/dy 



(6) 



and the condition d ^ IN^^'^ is simply Rc/d <C 1. It 
immediately hits the eye that HY's result is expressed in 
terms of purely classical quantities, given the particle's 
energy is equal to E^. Note also that the size of the ef- 
fect is different from KL's formula (^. For instance, it is 
much larger provided that Rc/d ^ W/Ej^. For a weak 
random potential, W <gC Ep ~ iJ^r, this inequality can 
be met simultaneously with Rc/d <^ 1. 

Naively, one might think that the weak levitation 
mechanism of HY stops functioning in lower magnetic 
fields where Rc ^ d. This turns out not to be true; 
however, the dependence of the weak levitation effect on 
magnetic field becomes slower. 



6En/E^ ^ {W/ENf{Rc/d) ^ {uct)- 



(7) 



[The final expression follows from r ^ {d/vF){EF /W)"^]. 
Equation (0) is the central result of this paper. It is rep- 
resented graphically in Fig. |l| together with the previous 
two. The plot should be understood as the dependence 
of the quantity SE^/E^ at the topmost Landau level 
on the ratio Rc/d. In other words, I assumed that the 
Fermi energy Ep is fixed but the magnetic field is chang- 
ing. For each value of the magnetic field one has to choose 
En closest to Ep. Of course, discreteness of N leads to 
some fine details on the curve in Fig. |l|. Such details are 
insignificant for large N, which is assumed to be the case 
for the most points on the plot. 

As one can see from Eq. (0), 5En /E'^ monotonically 
grows as the magnetic field decreases. Naturally, one 
would like to know if it ever becomes of the order of one. 
The answer is negative: even at the lowest magnetic field 
Bniin where the present approach is justified, the quan- 
tity SE^/E^ is still small. To verify that one needs to 
know what -Bmin is. Clearly, -Bmin is the largest of the two 
fields, at which the two simplifying considerations men- 
tioned in the beginning of the paper breaks down. One 



is the field Be where the separation into the cyclotron 
and drift motion ceases to be valid, and the other is the 
field where the quantum uncertainty (~ I) of the guiding 
center position becomes of the order of d. The crossover 
field Be was calculated in Ref. It corresponds to the 
point Rc/d ~ (Ep/W)"^^^ where the characteristic fre- 
quency of the drift motion becomes of the order of lUc- 
Combining this with the other condition, one obtains the 
largest value of the ratio Rc/d where the calculation is 
still valid, 



{Rc/d)r 



.{{Ep/Wf/'', kpd} 



(8) 



(kp is the Fermi wavevector in zero field). Substitut- 
ing this value into Eq. (Q) and keeping in mind that 
W <^ Ep, one obtains that ((5i?Ar/i?^)inax <■ 1 (see also 
Fig. 0). 
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FIG. 1. The sketch of the dependence of the relative 
size of the weak levitation effect SEn / at the topmost 
Landau level on the ratio Rc/d (see text). Labels "HY," 
"this work," and "KL" refer to Eqs. (g), (|), and @, 
respectively. Both axes are in log-scale. 

Let us now turn to the derivation of Eqs. and (|^). 
I am going to show that the effect is completely clas- 
sical and therefore the constraint W ^ fiujc imposed 
by HY is extraneous. The relevant condition is just 
Rc/d ^ {Rc/d)max- Also, in contrast to the quantum- 
mechanical treatment of Ref. |[ my instrument will be 
the classical perturbation theory. For this reason I will 
drop the subscript "|£." in Ep. To ensure continuity with 
the previous paper,E3 I will assume that the magnetic 
field is in the negative z-dircction, so that the cyclotron 
gyration is clockwise and the guiding center coordinates 
are given by 

Px = X + (Vy/Ujc), Py^y~ {Vx/Ujc), 

where x and y are the coordinate of the electron, and 
V = —V (sin 6*, cos 0) is its velocity. The most convenient 
form of the equation of motion for p is obtained if U is 
re-expressed in terms p and 9 only, i.e., as a new function 

V{p,0) = U[px + (v/ojc) cos d,py - (u/cJc) sin6']. 



where v — \/2{E — V)/m because of the energy conser- 
vation. The equation of motion for p is 
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dp/de = — -[zxVpVip,e)]. (9) 

This equation is of the Hamiltonian form with O/ujc, Py, 
mUcpx, and V being the time variable, the canonical co- 
ordinate, momentum, and the Hamiltonian function, re- 
spectively. The return to the original time variable can 
be accomplished by means of the equation 

= LJc + z[v X VU]/{mv'^). 

It is convenient to expand V{p, 9) in Fourier series, 

oo 

F(p,0)= (10) 

k— — oo 

If |fe| ^ Rc/d, the absolute value of Vfc is of the order of 
Wo = W(d/RcY^^; otherwise, it is much smaller. Sub- 
stituting Eq. ([l^) into Eq. (||), one obtains 

oo 

dp/d9= V [zxVpVkip)] (11) 

mujf. ^ — ' 

fe= — oo 

If one retains only the fc = term, then the right-hand 
side does will not depend on 9 and thus the guiding center 
motion will decouple from the cyclotron one. In this ap- 
proximation the guiding center performs the drift along 
the level lines of the potential Vq. In the d limit Vb 

is very close to the original potential U in agreement with 
the qualitative picture given above. If Rc is larger than d, 
then quite a few k terms are of the same magnitude 
as = one. In this case one can not simply ignore them; 
however, they can be made smaller by means of series of 
canonical (or almost canonical) transformations.E^ Each 
consecutive transformation reduces the oscillating terms 
by a factor of the order of 7 = Wq / {mui'^d'^) <C 1. In the 
end they become suppressed by a factor exp(— const/7). 
This program can be realized only if 7 1. Equation 
7 = 1 thus determines the magnetic field Be (see above) 
where the crossover from the adiabatic^rift to the ran- 
dom walk of the guiding center occurs £3 

For calculation of SEn to the first nonvanishing order 
in 7 only one such a transformation suffices. Let p and 
q be the new canonical coordinates after the transforma- 
tion. Define the "renormalized" guiding center coordi- 
nates, p^^ = pI (mujc) and p^'' = q. It is easy to see that 
p*^^^ has the following form 

P<')-p+^Ei[^x'^^^^] + 0(7'rf)- (12) 

The 0-independent term Vq in the Hamiltonian function 
is transformed into V^^ given by 

, 00 „ 

= + ^ V 4 i^^V-k X VVk] + 0(7'Wo). 

^ k=\ 

Define 



Uk{p.K) = j) ^e-'''^U[px + Rcoscb, Py + Rsincb], 

where K has the meaning of the kinetic energy and 
R = y^2K/muj^ of the corresponding cyclotron radius. 
Note a useful formula 

Uk{q,K) = i''e-''''''Jk{qR)U{q), (13) 

where tilde symbolizes the Fourier transform, q — 
q{cos 9k, sin 9k), and Jk is the Bessel function. 
It is easy to see that 

00 

(p, K)^Uo + J V(rfe - Zk) + Oij'Wo), (14) 

^ fe=i 

r. = |[vc/_.xvt/.], ^.^^ic/.r (15) 

The obtained expression agrees with the effective poten- 
tial of HY in the limit Rc ^ d. To see that one has to 
quantize the kinetic energy K = E — V°^ = fiujc{N +1/2) 
and keep only fc = 1 and k = 2 terms, which dominate 
the sum ia this limit. The levitation correction can be 
estimatedQ as 

SEr, ~ (y°*^)sp, (16) 

where "SP" stands for saddle-points, i.e., the points 
where 

(A) US = uy = 0, {B) us^uyy - ui^ur < o (i7) 

(the superscripts denote the partial derivatives) . The rest 
of the paper is devoted to the derivation of Eq. (0) from 
Eqs. ( p"4[]l6| ). I will assume that U(x,y) is an isotropic 
Gaussian random potential with zero mean. 

It turns out that for each k, Yfc and Zk are cor- 
related with at most one of the sets {U'q ,Uq} and 
{Uq^, Uq^ , Uq^}. Therefore, each time one needs to cal- 
culate either (Yfc) a or (Yk)B, or simply the unrestricted 
average (Yfc) (and similarly for Z^). The conditions 
and "B" are given by Eq. (|T7|). Notice that Yk and Zk 
are bilinear in U . This allows us to perform the (. . .) a 
averaging by means of the following general formula. Let 
Xi and X2 be linear in U , then 

{XiX2)a - {X1X2) - {usus)-^ 

X i{X,US){X2US) + {X,uy){X2Uy)) . (18) 

Similarly, 

{x,X2)b = {x,X2) + {us^un-'i{XiU'S''){X2US') 
i(XiC/-)(X2c/r) - Ux,uyy){X2usn)- m 

Equations (|8|) and ( p^ can be obtained using general 
properties of Gaussian potentials and isotropicity of U. 
The (. . .)a averaging is needed for Y2 and Zi, the (. . .)b 
averaging is needed for Yi, Y3, and Z2. All other Yfc's 
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and Zk's are to be averaged over the entire plane. The 
calculation is trivial but lengthy, and so I will give only 
the final result: 
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where Afc„ is defined as follows 



Ciq)MqR)JkiqR)q'', 



(20) 



with C being the correlator C(r) = (C/(0)C/(r)). If 
Rc » d, then 
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(21) 



This concludes the derivation because Eqs. (^^ and (|21 
immediately give Eq. (0) . 

The correction (y° — Vq) has quite peculiar proper- 
ties: it is typically positive at the points where the gra- 
dient squared (Vt/p)^ is smaller than its average value, 
typically negative otherwise, and almost vanishes on 
average. Indeed, one can show that the unrestricted 
spatial average of is {V^) = -Aoi/{2mujlR) - 
-{V°^)sY>{d/Rcf, which is much smaller than {V''^)sy> 
by the absolute value. 

I speculate that the last property becomes important 
in the vicinity of Be where the crossover from the drift 
to the diffusion occurs. In the diffusive regime the tra- 
jectory of the guiding center is no longer bound to the 
level line V^^ = const but samples the entire area. Thus 
(F°'^)gp should approach {V^^). The latter is indistin- 
guishable from zero within the accuracy of such an un- 
rigorous argument. Thus, I expect the ultimate downfall 
of the solid curve in Fig. ^ near its termination point. 

Concluding this paper, let us emphasize that the dis- 
crepancy between KL's formula (^) and Eqs. (^) and (|^) 
does not contradict to Eq. (Q), which comes from the 
scaling arguments. Indeed, for B < Be the "classical" 
or the "unrenormalized" Hall conductance is determined 
not by the average density n but by the density Up near 
the percolation contour which is the area responsible for 
the transport; therefore, 



(22) 



' xy — '^p- 
mUJc 



Equation ( p^ now follows from Eqs. (|^), (p2[), and 
Up = [m/Tih ){Ep — V°^) ("p" again means "percola- 
tion"). At B > Be the percolation contour is of no im- 
portance so that Eq. (^) becomes valid, and presumably 
so does Eq. (^ as well. 

In conclusion, I showed that the relative size SEn /E"^ 
of the weak levitation effect is always much smaller than 
one. Moreover, it is expected to decrease near its termi- 
nation point B — Be- This strongly suggests that this 
effect has nothing to do with the strcmg levitation pre- 
dicted by KhmelnitskitI and Laughlin.BI 



Finally, I should also mention that SEn can be mea- 
sured experimentally by chartingi-the global phase dia- 
gram of the quantum Ifall effect, til the procedure initi- 
ated by Glozman et al.tB Very interestip|g-|and puzzling 
findings have been reported recently.ll3il3 However, it 
seems that the electron-electron interaction plays a cru- 
cial role in the observed phenomena. This complicated 
issue is beyond the scope of the present paper. 
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